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We analyze the general conditions on the equation of state w required for quantum fluctuations of 
a scalar field to produce a scale-invariant spectrum of density perturbations, including models which 
(in the four dimensional effective description) bounce from a contracting to an expanding phase. 
We show that there are only two robust cases: ui ~ — 1 (inflation) and lo S> 1 (the ekpyrotic/cyclic 
scenario). All other cases, including the w ~ case considered by some authors, require extreme 
fine-tuning of initial conditions and/or the effective potential. For the ekpyrotic/cyclic {w ^ 1) 
case, we also analyze the small deviations from scale invariance. 



Until recently, the only known mechanism for generat- 
ing super-horizon, adiabatic, scale-invariant density per- 
turbations was inflation [1] . One of the great surprises of 
the ekpyrotic [2-4] and cyclic models [6,7] is that they 
provide a second, alternative mechanism. In this pa- 
per, we introduce a gauge invariant, systematic analysis 
(including gravitational backreaction) that identifies the 
most general conditions required assuming the density 
perturbations arise as a result of quantum fluctuations 
of a single scalar field with potential V{(j)). We in- 
clude the possibility that the universe may bounce from 
a contracting to an expanding phase, and that the per- 
turbations can be matched across such a bounce in an 
unambiguous way [3-5,8]. 

First, we must obtain scale-invariant fluctuations dur- 
ing a period in which the scalar field dominates the en- 
ergy density of the universe. The conditions for this to 
occur can be characterized by the equation of state w dur- 
ing this epoch. We find three interesting cases, for each 
of which w is nearly constant: i) an expanding universe 
with w ~ —1, corresponding to slow-roll inflation; ii) a 
contracting universe with w 3> 1, corresponding to the 
ekpyrotic/cyclic models; and, in) a contracting universe 
with w « 0, as discussed by Wands [9] and by Finelli and 
Brandenberger [10]. Although the last case does gener- 
ate a scale invariant spectrum of curvature perturbations, 
we shall show that the corresponding Newtonian poten- 
tial has a very red power spectrum. This points to a 
serious instability of the background, which we explic- 
itly identify in the infinite wavelength limit. Unlike the 
ekpyrotic/cyclic cases, the w — background solution 
is unstable and hence not a dynamical attractor. The 
only additional cases are ones in which w is rapidly time- 
varying, but these require extreme fine-tuning of V{(l)). 

Hence, we find that the ekpyrotic/cyclic case remains 
as the only viable alternative to inflation, assuming suc- 
cessful matching of the growing mode perturbations gen- 
erated during the contracting phase onto growing mode 
perturbations in the expanding phase, as proposed in 
Refs. 4, 5, and 8. Using this prescription, we derive 
an expression for the deviation from scale-invariance for 
cyclic/ekpyrotic models in terms of "fast- roll" parame- 



ters e and ry, analogous to the "slow-roll" parameters of 
inflation. 

Our analysis is restricted to the case of a single scalar 
field, which includes the simplest inflationary scenarios as 
well as the ekpyrotic/cyclic models. In Newtonian gauge, 
the perturbed metric for a spatially-flat background can 
be expressed in terms of a single gauge invariant variable 
$, the Newtonian potential, as 

ds^ = a2(r) •{-(! + 2$(f,T))dT2 + (l- 2$(f,r))rff2}, 

where r is conformal time, and where we have used the 
fact that fluctuations of a scalar field do not generate 
anisotropic stress. 

While knowledge of $ is sufficient to determine the 
perturbed metric, it is useful to introduce a second vari- 
able, C, which is the curvature perturbation on comoving 
hypersurfaces [11,12]. $ and C are related by 



^ 3a^{l + w) \a'/a^ J ' ^ ^ 

where a prime denotes differentiation with respect to r. 
The variable C has the virtue of remaining nearly con- 
stant at snperhorizon wavelengths during epochs of ex- 
pansion. In inflationary models, it allows one to eas- 
ily match the Newtonian potential at horizon reentry in 
the matter dominated phase to that calculated at hori- 
zon exit during the inflationary phase. In the ekpy- 
rotic/cyclic scenarios, once its spectrum has been deter- 
mined after the bounce to an expanding phase, ^ also gives 
the perturbation amplitude at horizon re-entry. 

It is tempting to suppose that we should only be in- 
terested in tracking the evolution of ^. However, even 
though ^ is continuous throughout a contracting phase 
or an expanding phase, it can undergo a rapid jump dur- 
ing the transition between the two [4,8]. It is necessary 
to match the incoming <I>i„ and Cj„ to the outgoing ^out 
and Cotit- Generically there is some mixing and Cout de- 
pends on a combination of the <I>i„ and [4]. Hence, it 
is important to know both and Qn at the bounce. 

For the ekpyrotic/cyclic case, $i„ is scale- invariant and 
Qn is blue (decreasing at long wavelengths), and, hence. 
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Cout is dominated by the scale- invariant contribution (due 
to at long wavelengths, leading to a scale- invariant 
spectrum as modes re-enter the horizon during the ex- 
panding phase. For the w ss case, Qn is scale-invariant 
before the bounce and $j„ turns out to be red (increasing 
at long wavelengths and exponentially larger than Qn). 
This suggests that the background possesses a serious 
long wavelength instability; we confirm this by showing 
that the background solution is not a dynamical attrac- 
tor. 

The $ spectrum: The differential equation for the k- 
mode Ufc of the gauge invariant variable u, related to the 
Newtonian potential by m = a^/<j)' (where we drop the 
(in) subscript, henceforth), is 



uk = 0; 



P{T)=T^H^aUe + 




(3) 



where H = a'/a^ is the Hubble parameter, N = In a, 
and e = 3(1 -|- w)/2. Note that in the case of inflation, 
e reduces to the usual "slow- roll" parameter, while N is 
the number of e-folds of expansion. 

In the regime fc^r^ » \(3\, Eq. (2) reduces to the equa- 
tion for a simple harmonic oscillator, and Uk is stable. 
When k'^T^ <C \/3\, however, the amplitude of the mode 
is unstable. In order to have a situation where successive 
modes with increasing k are becoming unstable and grow- 
ing, we need (3 positive, and assuming that (3 is slowly 
varying, we require that r be negative and increasing. 
This applies to expanding models, such as inflation, or 
contracting models, such as the ckpyrotic/cyclic scenar- 
ios. 

For general time-varying w, (3{t) will be a complicated 
function of time, and one can use numerical methods to 
solve Eq. (2). In the most plausible cases, however, it is 
reasonable to approximate w as constant, at least for the 
observationally relevant range of modes. 

It is well-known [4,13] that solutions with constant w 
correspond to potentials of the exponential form, = 
—Voe~^'^, where c and Vq are constants. In this case, the 
equation of state is related to the slope of the potential 
by 



e= -il + w) 



\ Ya 

2\V 



(4) 



(we use units where SttG = 1) and the solution for the 
background, assuming homogeneity, isotropy and spatial 
flatness, is given by 



a(T) ~ (-t)^/(^"-^\ H 



1 



e- 1 



V{<t>) 



l)ar 



3) 



(e- I)2a2r2' 



(5) 



Substituting the above into Eq. (3), we obtain (3 = 

e/(e — 1)2 . Since /3 is constant in this case, Eq. (2) can 
be solved analytically, with general solution 



Wfe = V^[Ci(fc)J„(-fcT)+C2(fc)J_„(-fcT)], (6) 



where n = \/ P -\- 1/4, J„ is the Bessel function of the 
first kind of order n, and Ci{k),i = 1,2 are arbitrary 
functions of k. 

The functions Ci{k) are determined by specifying ini- 
tial conditions when the mode is stable, i.e., when 
(2) fc^T^ » (3. In this limit, we make the usual assumption 
that the fluctuations in (f) are in their Minkowski vac- 



uum, which corresponds to Uk 



^^/(2A:)3/2. Using 



the relation Uk = a^k/fj^' and the asymptotic properties 
of Bessel functions, this gives 



Uk 



\/—TTkT 

4A;^/2 sin(7rn) 



[J_„(-ftr)-e-'""J„(-fcT)], (7) 



where we have neglected an irrelevant phase factor. 

We are interested in the amplitude of the mode in the 
long-wavelength regime, fc^r^ <^ f3. In this limit, we can 
expand the Bessel functions to obtain 



23/2 sin(7rn)r(l - n) 

2n 

1(1 — ri ] I — ft"r \ 

1 - e"""* 



-n-l-l/2 



,r(l -n) ( -kT 



r(i + n) 



r(l -n) ( -kT 
r(2 - n) 



(8) 



The Newtonian potential has a scale-invariant spectrum 
if the rms amplitude of ^k varies with k as /c~3/2_ Hence, 
we conclude that this will be the case if n « 1/2. Re- 
calling that n = \/ /3 + 1/4, this can be expressed as a 
constraint on /3 or, equivalently, e: 



< 1. 



(9) 



Wc have thus translated the requirement of scale in- 
variance for $ into a condition on the background equa- 
tion of state. Therefore, we may now determine what 
choice of w will satisfy Eq. (9) and lead to a Harrison- 
Zel'dovich spectrum. 

First, this condition is clearly satisfied when e <C 1, 
that is, when w « — 1. This corresponds to the case of 
slow- roll inflation [1]. Note, however, that there is a sec- 
ond regime in which condition (9) holds, namely when 
e ^ 1, corresponding to w 3> 1. This is the limit rel- 
evant to the production of fluctuations in the ekpyrotic 
and cyclic scenarios [2-4,6,7]. These two regimes are in 
some sense at opposite ends of parameter space. In the 
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inflationary case, e plays the role of a slow-roll parameter 
and is therefore small. In the ekpyrotic and cyclic sce- 
narios, however, e is large compared to unity. Also, from 
Eq. (5) ; wc sec that the universe is expanding in the first 
case and contracting in the second. 

This analysis assumed a nearly constant w so that /3(r) 
in Eq. (3) and, consequently, the spectral index is nearly 
constant. Note that this assumption is not necessary. It 
is possible, in principle, to build models for which the 
time-variation of w is non-negligible, and yet the deriva- 
tive terms in Eq. (3) conspire to cancel for a significant 
range of e-folds, N. This has been discussed for inflation 
by Wang et al. [14] who showed that maintaining the can- 
cellation for many e-folds requires extreme fine-tuning of 
V{(p) compared to the constant w cases. Hence, these 
models seem highly unlikely. 

The C spectrum: To calculate the spectrum of the sec- 
ond gauge invariant variable of interest, we substitute 
the expression for $ obtained in Eq. (8) into Eq. (1). The 
leading term in the expansion for $ is 

(10) 

where we have omitted the numerical coefficient. As be- 
fore, we approximate w as constant, and, thus, so is n. 
Using Eqs. (5), we find 

as well as 

^ ^ (_^)-(l+e-)/(e--l) . (12) 

Since the expression for A involves a factor of |e— 1|, 
we must consider the cases e > 1 and e < 1 separately. 
Starting with the latter (which includes slow-roll infla- 
tion), we find A = 0, and therefore 

fc3/2 ( J^") ^ fc-.7(l-e-)(_^)(l+?)/(^-l) . (13) 

\a' / a-^ J 

Using the relation between ( and $ given in Eq. (1), we 
see that ( will have a scale-invariant spectrum if e <C 1. 

Recall that this limit corresponds to slow-roll inflation 
and that it also leads to a scale-invariant spectrum for 

For the case e > 1, which includes the ekpyrotic/cyclic 
scenarios, we find A = — (l-|-e)/(e— 1), and thus 

A:3/2 ( ^ ^-i/(,-i) _ (14) 

Since the right hand side is independent of time, Eq. (1) 
implies that this leading term for $ does not contribute 



to C. In order to determine the long-wavelength piece of 
^, we must therefore keep the higher-order terms in the 
expansion for $ given in Eq. (8). It is straightforward to 
show that the result is of the form 

k^^\ ~ fl{T)k'/^'-^^ + /2(T)fc(2^"-3)/(^"-l) , (15) 

where /i (r) and /2 (t) are time-dependent factors. 

For the ekpyrotic/cyclic scenarios, corresponding to 
the regime e :» 1, the first-term in Eq. (15) gives the 
dominant contribution at long- wavelengths, and thus 
^3/2^ goes like k. Hence, while the condition e 1 led 
to a scale-invariant spectrum for $ in the pre-big bang 
phase, it yields a blue spectrum for (. As we have seen 
above, this is a consequence of the fact that the growing 
mode of which is scale-invariant in this limit, is pro- 
jected out of (. Thus, C is determined by the next-order 
correction in the expansion for $, which is down by a 
factor of fc^" « k. Nevertheless, as mentioned earlier, 
C jumps at the bounce (mixes with the scale-invariant 
$ mode) and is dominated by the scale-invariant contri- 
bution at long wavelengths after the universe begins to 
expand. 

Note that, if we choose e w 3/2, then the second 
term in Eq. (15) dominates at large wavelengths, and 
the resulting spectrum for ( is nearly scale-invariant be- 
fore the bounce. This case, identified by Wands [9] and 
recently studied by Finelli and Brandenberger [10], de- 
scribes a contracting universe with a dust-like equation 
of state, w w 0. However, using our results above, we 
see that the corresponding spectrum of $ is strongly red 
(fc3/2^^ ~ k~'^T~^), indicating a severe long wavelength 
instability. Then if $ and ^ mix at all at the bounce, 
the red contribution would dominate at long wavelengths 
rendering the resulting universe phenomenologically un- 
acceptable. 

Stability of the background solution: Before consider- 
ing evolution of quantum fluctuations, we should first 
consider whether the constant w solutions we have as- 
sumed as background solutions are stable attractors to 
the equations of motion. If they are not, then they could 
not arise in a cosmological solution without extraordinary 
fine-tuning of initial conditions. The expanding inflation- 
ary {w = —1) is known to be a stable attractor. Here 
we show that the contracting ekpyrotic/cyclic (th ^ 1) 
phase is also a stable attractor, but the contracting w = 
phase is not. 

The stability of the background solution may be stud- 
ied in the inflnite wavelength (i.e., homogeneous) limit 
simply by considering the scalar field equation in a ho- 
mogeneous Universe: 

^ + 3-<^ = -y^, (16) 
a 

with dots denoting derivative with respect to proper time 
t. The Priedmann equation allows us to express a/ a in 
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terms of the scalar field. (For simplicity wc do not per- 
turb the space curvature since the above calculation in- 
dicates this effect is subdominant.) Setting = (j)B+S(j), 
with (j)B the background "scaling solution," 

a{t) oc i-tr, V{cI>b) = -Voe-^^^W = ^(3^ _ iyt\ (17) 



where p = (16) becomes 



S4> + 



t 



l-3p 



= 0, 



(18) 



whose two linearly independent solutions for p ^ | are 

5(f> ~ and t^^^^. For all the first solution is just an 
infinitesimal shift in the time to the big crunch: S(f> oc 4>b- 
Such a shift provides a solution to the Einstein-scalar 
equations because they arc time translation invariant, 
but it is physically irrelevant since it can be removed by 
a redefinition of time. In contrast, however, the second 
sohition represents a physical perturbation of the back- 
ground solution. For p > 5 , the second solution grows as 
t approaches zero, indicating an instability of the back- 
ground scaling solution. (For p = |, the second solution 
is ln{—t).) This is confirmed by calculating the ratio 
of kinetic energy to potential energy in the scalar field, 
which is constant and equal to 



K 
V 



1 



w 



1 



1 



3p- 1 



(19) 



in the background solution. This ratio is unaltered by 
the first perturbation solution (since it is a time shift) 
but the second yields 



3p 



(20) 



which for p > ^ diverges as t = is approached. This 
shows in particular that the w; = (p = |) background 
scaling solution is unstable and hence not an attractor. 
Conversely, the ekpyrotic/cyclic cases, which correspond 
to p << 1, possess scaling solutions which arc are sta- 
ble attractors in the infinite wavelength limit, since the 
only growing mode is as we have discussed just a time 
translation. 

A more general, if more heuristic argument can be 
obtained by comparing the Priedmann equation for the 
three cases, 



1 



(21) 



where H is the Hubble parameter and p is the energy 
density. For the expanding w = —1 case, the scalar field 
energy density is nearly constant, whereas the curvature, 
radiation (oc 1/a**), matter (a 1/a^), and other forms of 
energy density all decrease as a expands. Hence, the con- 
tant energy density w = —1 state is an attractor. For the 



contracting w ^ 1 case (a ^ t^ with p ^ I), the scalar 
field energy density increases as Xja}!^ as a decreases, 
whereas the curvature, matter, radiation, or other forms 
of energy density increase at a slower rate. Hence, the 
w ^ \ contracting solution is also an attractor of the 
Priedmann equation. However, for the w w solution, 
the scalar field energy density increases as 1/a^, but ra- 
diation density increases at a more rapid rate. Further- 
more, as shown above, a small perturbation drives the 
universe away from w « (K « V^) towards w = \ 
(K y), a state in which the energy density (ex 1/a^) 
increases even more rapidly. Hence, the w ^ back- 
ground is not a stable attractor. 

We can state the conclusion more generally. For 
an exponential potential proportional to exp(— C(?i), we 
have a scaling solution when V is positive (negative) for 
c < (c > a/H). In the scaling solution, we have 
w = '^'/{'ip) — 1 = (c^/3) — 1. It follows from the analysis 
above that in a contracting universe, the scaling solu- 
tion is only a stable attractor if w > 1 (or c > \f%) and 
the scalar potential is negative. By time reversal, we in- 
fer that in an expanding Universe the scaling solution is 
only a stable attractor if w < 1 (c < ^) and the scalar 
potential is positive. Inflationary and 'quintessence'-type 
scaling solutions are both included in this latter case. 

Thus, we have completed our classification of the pos- 
sibilities. In particular, we have seen that when the uni- 
verse is contracting with w ^ 1, as in the ckpyrotic 
and cyclic scenarios, the Newtonian potential $ devel- 
ops a scale-invariant spectrum while that of C is blue. 
Provided <I> and C. mix at the bounce (which is to say 
that the growing mode of the contracting phase does not 
match to a pure decaying mode in the expanding phase), 
one obtains scale-invariant density perturbations in the 
expanding phase. If any radiation is generated at the 
bounce, mixing is expected [4,8]. When the universe is 
contracting with a dust-like equation of state (w « 0), C 
acquires a scale-invariant spectrum, while <f> acquires a 
red spectrum. With mixing at the bounce, one obtains 
an unacceptable red spectrum. More generally, we have 
shown that the w — scaling solution is not an attrac- 
tor. Therefore, there is no reason to expect the universe 
to reach the scaling solution in the first place. More gen- 
erally, we have shown that the only way to obtain a sta- 
ble attractor scaling background solution in a contracting 
universe is to have a negative scalar field potential, as in 
the ckpyrotic and cyclic models. 

Spectral index in ckpyrotic and cyclic models: In the 
remaining part of this paper, we shall focus on the ekpy- 
rotic/cyclic generation of perturbations and calculate the 
spectral index, giving a treatment analogous to that given 
for inflation by Wang et al. [14]. 

Recall that approximate scale invariance of the power 
spectrum in the ekpyrotic/cyclic scenario requires that 
e be large and nearly constant as modes become unsta- 
ble. Since e 1, it is convenient to introduce a small. 
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"fast-roll" parameter e as 



(22) 



The condition e ^ 1 implies e ^ 1, which translates into 
the requirement that the potential be steep. 

Since e is large and nearly constant, the parameter /? 
defined in Eq. (3) reduces to 



1 /dine 



2\ dN 



(23) 



where we have assumed that rf^ Ine/dN^ and {dlne/dN)^ 
are much smaller than d\ne/dN. 

Recalling from Eq. (4) that e Ri V^/2V^ for nearly 
constant e, we obtain 



rfln 



d\n i 



dN \aHj # -H^Jl^J^' ^''^ 

where we have introduced a second fast-roll parameter r], 
defined by 



^- VI 



(25) 



Note that = corresponds to pure exponential poten- 
tials. 

Substituting for (j)' /{aH) using Eqs. (5), Eq. (24) re- 
duces to 



dine 
dN 



4er). 



(26) 



Since e is assumed to be nearly constant and large, 
Eq. (26) implies \r]\ <^ 1; that is, the potential must be 
nearly exponential. 

Prom the background solution given in Eqs. (5), it is 
easily seen that 



THa^i^ •{l + 0(e,r?)}. 



(27) 



Therefore, substituting Eqs. (22), (26) and (27) into 
Eq. (23), we find 



/3«2(e + 7?). 



(28) 



We may now proceed to calculate the spectral index 
of density perturbations. As seen from Eq. (8), the long- 
wavelength limit of the Newtonian potential is given by 



(29) 



where we have used the fact that n = P + 1/4 and 
/? ^ 1 in this case. 

Equation (29) describes the spectrum of ^ for modes 
that went unstable during the contracting phase. In 
Ref. [4], it was argued that the pre-big bang spectrum 



of $ gets imprinted on the long-wavelength part of C as 
the universe undergoes reversal from contraction to ex- 
pansion. (See also Refs. [7], [8] and [15].) Then, the 

post-big bang spcctnim of energy density perturbations 
is given by 6k ~ k^^ , corresponding to a spectral index 



1 



dln\Sk\' 
dink 



-2/3, 



(30) 



where rig = 1 corresponds to an exactly scale-invariant 
(Harrison-Zcl'dovich) spectrum. 

Substituting Eq. (28) into Eq. (30), we find that the 
spectral index of density perturbations in the cyclic and 
ekpyrotic scenarios is given by 



n. - 1 



-4(e-Fr?) 




In the case of pure exponential potentials, r] vanishes 

identically, and therefore the spectrum is red (since f > 
0). For potentials of larger curvature than an exponen- 
tial, such as —e~'"^ , one has ?7 > and the spectrum 
is also red. However, for potentials of smaller curvature 

than an exponential, such as e~'^V^, one has ry < 0, and 
the spectrum will be blue if e+r] is also less than zero. For 
instance, the string-inspired potential of Ref. [2] led to a 
blue spectrum. We now see that both red and blue spec- 
tra can be achieved, as anticipated by Linde et al. [16]. 

It is instruc;tive to compare Eq. (31) with its counter- 
part in slow-roll inflation [14,17] 



1 = -6e + 2fj, 



(32) 



where e = V'^/2V'^ and fj — V,p,p/V arc the usual slow- 
roll parameters of inflation. It is easily seen that pure 
exponential potentials also yield a red spectrum in this 
case. Once again, it is possible to find potentials for 
which the spectrum can be either red or blue. 

In summary, our work shows that the inflationary and 
recently introduced ekpyrotic mechanisms arc the com- 
plete set of approaches for generating a nearly adiabatic, 
scale-invariant spectrum of fluctuations from a single 
scale field without extreme fine-tuning. For the ekpyrotic 
case, we have shown how the spectral index is related to 
fast-roll parameters that characterize the slope and cur- 
vature of the scalar field potential scale. In Ref. [7], we 
show that this constraint requires essentially the same 
amount of fine tuning as the slow-roll conditions for in- 
flation. In Ref. [18], we consider the spectrum in a mixed 
case where the scalar field rolls from an expanding infia- 
tionary regime to a contracting ekpyrotic regime. 
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